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Abstract 

A simple mechanism for the confinement of red cells in the middle of narrow blood vessels 
is proposed. In the presence of a quadratic shear, red cells deform in such a way to loose fore- 
aft symmetry and to achieve a fixed orientation with respect to the flow. This leads to a drift 
away from the vessel walls, when the vessel diameter goes below a critical value depending on the 
viscoelastic properties and the dimensions of the cell. 



PACS numbers: 87.45.Hw, 47.15. Gf, 47.55.Kf, 



Blood is essentially a concentrated suspension of red cells in plasma. Other particles, like white 
blood cells and platelets make up less than 1% of the particulate. Due to the high volume fraction, 
(around 40% — 45%) of the red cells, and to their mechanical properties, blood rheology exhibits 
rather complex behaviors. Focusing on steady state situations, in particular, blood viscosity turns 
out to depend on the diameter of the vessel in which it is flowing. 

In large blood vessels, with diameters above ~ 100/ini, particles are randomly distributed 
in the plasma; the effective viscosity appears to be constant, and could be calculated, to a first 
approximation, using Einstein theory 0|. For diameters between lO/im and 100/im, red cells align in 
the middle of the vessel, resulting in a decrease of the effective blood viscosity, which is fundamental 
for the well functioning of the cardiovascular system: the so called Fahraeus-Lindqwist effect . 
For still smaller diameters, corresponding to the capillary range, these cells have to squeeze their 
way through the vessel, leading to a sharp increase in the effective viscosity. In all three cases, red 
cell deformability plays an important role, facilitating the particle flow and, as it will be shown 
below, contributing to the alignment process in the interediate range. 

The alignment of red cells in the middle of small blood vessels, corresponds to the intuition, that 
particles in suspension should strive to minimize dissipation and to accommodate themselves in 
the regions of minimum shear. A dynamical description, however, is not simple and requires either 
consideration of inertia ^ , or of the cells non-sphericity |^ and deformability, and possibly, of 
their mutual interactions. 

The purpose of this letter is to investigate the hydrodynamic interactions that stay at the 
basis of the Fahraeus-Lindqwist effect and to propose a simplified model for the description of the 
process, based mostly on geometrical reasoning. It appears that the drift produced by inertial 
corrections (see e.g. |6)) cannot account for this effect. Deformation by the quadratic part of the 
shear, turns out instead to be the essential ingredient for the migration of a red cell in a small blood 
vessel. This confirms results of numerical simulations carried out in |Q, and, more recently, in 
. Interaction among the red cells contributes to the effect but is not essential to the migration 
process. 

Consider a particle in a channel flow, which, in a reference system with origin at the particle 
center, can be written in the form: 

V = e3(a -I- /3a;2 -f- 72^2)- (1) 

It should be mentioned that channel flow configurations are not just a device for the derivation of 
simpler models, but have physiological relevance in the description of blood flows in internal organs 
such as the lungs. The presence of the particle produces a perturbation in the velocity field, which 
is determined imposing no-slip boundary conditions at the particle surface and at the channel 
walls. If the channel gap and the particle are sufficiently small, the fluid will be described to first 
approximation by Stokes equation. Thanks also to linearity of this equation, a qualitative idea 
of the velocity perturbation can be obtained, writing it as a superposition of images and counter- 
images analogous to those of electrostatics: v-|-v' -l-v" -I- ... Here, v is the velocity perturbation in 
the infinite domain case, is the first image, obtained enforcing no-slip on the walls: v -|- v' = 0, 
v" is obtained imposing no-slip on the particle: v' -I- v" — 0, and higher orders are obtained 
similarly, enforcing no-slip alternatively on the walls and on the particle. 

If the particle is fore-aft symmetric, expanding the velocity in vector spherical harmonics JTo[ , 
it is easy to see that v~ ~ a/?i?(i?/Z)^^ -I- b^B?{R/lY'^~^^ , where R is the particle size and / is the 
distance to the closest wall. If the ratio R/l is sufficently small, the drift of the particle is given 
approximately by the value of at the particle center. This allows to visualize the mechanism for 
lift close to a wall, as shown in Fig. I. 

It is well known that spherical symmetry does not allow for transversal drift in creeping flow 
conditions. Inclusion of inertia is sufficient to break the symmetry of the velocity lines (see Fig. 
fa), and this leads to a transverse migration of the order of Rep(3R with Rep = the particle 
Reynolds number and 1/ the kinematic viscosity of the fiuid. Due to the smallness of Rcp, however. 
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Figure 1: Sketch of velocity perturbation around a particle in various shear flows, and image 
velocity field produced by a wall; both velocity fields are in the particle's reference frame, so that 
the wall is moving upwards, (a) effect of inertia: the velocity lines have an increasingly stronger 
tilt along the flow, as one gets farther away from the particle; (b) perturbation by an ellipsoidal 
particle in linear shear, in the absence of inertia; (c) the same in a quadratic shear: the velocity 
V has an octupole rather than a quadrupole symmetry, but the same drift of the previous case 
occurs. 



inertial drifts turn out to be too small to produce a Fahraeus-Linqwist effect. 

A second mechanism for symmetry breaking lies in the non-sphericity of the cells. The drift 
of ellipsoidal particles has been studied in the context of sedimentation More recently in 

following in part jl^, the effect of tank-treading motions in linear shear has been taken 
into consideration. The velocity perturbation produced by an ellipsoid, is not symmetric even to 
0{Rep), and a transverse drift results, as illustrated in Fig. Ib-c. In the case of a linear shear, the 
velocity perturbation v decays quadratically at large distances . This leads to the expression for 
the drift in the presence of a wall ||, p^ : 

v\l)~C{9,R/l)l3Rsm29, (2) 

where C{e, R/l) = CqR^ + cR^/l* + c-yR-^/ipP) + with Co independent of 0. The 0{{R/lY) 
part of the correction is basically due to the near field of v, while the 0{R'^ / {(il^)) contribution is 
due to the quadratic part of the shear. Analysis carried out in shows that the lift in a channel, 
for R/l small, is obtained to a very good degree of confidence, by linear superposition of the drift 
produced individually by each of the walls. Thus, particles will be drifting at maximum velocity, 
away from the wall, when the longest particle axis is in along the stretching direction of u, and 
towards the wall, when it is along the compressing direction. 

A tank-treading cell could maintain a fixed orientation, thus leading to a steady drift. Unfortu- 
nately, due to their stiffness, red cells will resist tank-treading motion in vivo iQ, and will rotate 
as rigid objects, under the effect ot the vorticity part of the shear ("flipping" motion regime). 

In order for migration to take place, it is necessary that either the particle undergoes a kind 
of rigid rotation which privileges orientations corresponding to drift away rather than towards the 
wall, or that some new mechanism for fixed orientation is present. At rest, red cells can be roughly 
described as axisymmetric oblate ellipsoids, with major axis ~ 8/xm and minor axis ~ 2yum [ p!5| . 
Introducing a reference system {x'lX^x'^}^ with x'l = xi, and x'^,x'2 respectively along the minor 
and the major ellipsoid axes, the external shear can be written: 

e3a;2 = ;^ (632:^2 - ^'2^'^) + + e'jXg) cos26' -f ^(62X3 - 63X3) sin 26* (3) 



2 

GaXj = cos 0[e3(cos^ 9x2^ + arc? 9x'^^) — 63 sin 2^X2X3] 



and 
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+ sin6l[e2(cos2 6*4^ + sin^ - e'^ sin 2ex2x'3] (4) 

Because of symmetry under reflection across the planes x'iX2 and 2^1 Sg, the term in sin 20 in Eqn. 
(3) does not contribute to torque. Turning to the quadratic shear, all terms entering Eqn. (4) are 
symmetric either under reflection: x[ — > — x^, with i — 2,3, or simultaneously: Xj, X3 — > —x'2, —x'^- 
Hence the quadratic part of the shear does not produce any torque on an ellipsoidal cell and does 
not contribute to the angular velocity 9. The contribution to 9 produced by the linear part on a 
fore-aft symmetric cell, instead, has the form pc| ]: 

9 = P{a-bcos2e), (5) 

Hence, the cell will rotate with an even angular velocity: 9{6) = 9{—9), while, from Eqn. (2), its 
center will undergo a zero mean oscillatory motion in the X2 direction. For an isolated cell, it is 
then necessary to take into consideration the symmetry breaking deformations produced by the 
external flow. 

The linear shear causes a modification of the cell eccentricity, which leads to corrections pro- 
portional to {AR/ R) sm29, with AR the deformation, in the coefficients b and Co [see Eqns. (2) 
and (5)]. For a membrane shear elasticity k, one has dimensionally: AR/R ~ pv(3R^ . 

The quadratic shear in Eqn. (4) is in the form ucos6' + wsin6'. Considering a fixed cell, the 
two terms u and w produce respectively deformations which are even in x'2 and x'^^ but break 
fore-aft symmetry in x'^ and x'2. These deformation will couple with w and u to produce a torque. 
Hence, for small deformations, this torque will be in the form: Mi = sin26'/(6'), with f{9) bilinear 
in u and w. From dimensional analysis, for the kind of flow described by Eqn. (1) and small cell 
deformations, one finds 

e^DP^,u.29, (6) 

where p and k are the plasma density and the membrane shear elasticity. This torque has the 
odd dependence on 9, necessary to break the symmetry between cell orientations associated with 
inward and outward drifts. The dimensionless coefficient D depends on the cell rest shape and on 
the angle 9, with this last dependence disappearing when the rest shape is almost spherical [p^ . 
Notice that the key ingredient for the production of this torque is the violation of fore-aft symmetry. 
For this reason there is no coupling between quadratic shear and the deformation produced by the 
linear strain, and also the viceversa is true 

A tractable model can be obtained if one supposes that the elastic time scale is faster than that 
for fluid motion, which is equivalent to assuming small deformations. (Smallness of the deformation 
is necessary also to make sure that interaction with the walls are the only source of drifts in the 
problem [|o|). 

Consider then a channel flow of width 2d and fluid velocity at the center V (see Fig. 2). After 
expressing all lengths and times in units of d and d/V , the model equations will read: 



2r^x{C{9, l-x)- C{9, 1 + x)) sin 29 
2x{a-b cos 261) + Dr^Vpi^n-'^ sin 29 



where r ^ R/d and x — X2/d are the dimensionless particle radius and distance from the vessel 
axis. Following Eqn. (2), the two terms C in the first line of Eqn. (7) account for the transverse 
drift, as sum of the individual contributions from the two walls. The factor x in front gives the 
strength of the linear part of the shear at the particle center. The second line contains the balance 
between the torque from the vorticity part of v (even in 9) and the counter-torque (odd in 9) 
provided by the interaction between quadratic shear and the cell deformation. 

Taking D independent of 9, as in the spheroidal limit, allows to analyze the dynamics of the 
system (actually, it is enough that D does not change sign). Its phase portrait is shown in Fig. 
3 and clearly indicates the presence of a stable fixed point a,t x — 0,9 = + n)TT corresponding 
to red cell alignment in the channel axis, perpendicular to the flow. However, if r is small, the 
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Figure 2: Ellipsoidal cell in a channel flow. At the particular position and orientation represented, 
the cell is moving downwards while rotating clockwise (i.e. 9 < 0). 



approach to the channel axis will be very slow, as it is easy to see from the following considerations. 
The first of Eqn. (8), for r small and away from the wall, can be approximated as i ~ r^x^ sin 20; 
substituting into the second one and expanding in powers of r one finds: 

din a; sin 20 / Dr^ sin 29 \ 

~d9 a-6cos20V"^~ 2a;(a - &cos20) ~^"') ^' 

and the drift will be 0{r^). Terms of the same order are produced including the corrections to b and 
Co produced by the deformations from the linear shear. Thus, the fixed point basin of attraction 
has an effective boundary at a; ~ r^VpiyK~^ , the width of the region of x occupied by the curve of 
turning points for (the dotted line in Fig. 3). An "isolated particle" Fahraeus-Lindqwist effect 
will be present only when the width of this region and the one of the channel are comparable, and 
this leads to the condition on d: 

d<E(^)^R, (9) 

with E depending on the rest shape of the cell through the coefficients a, 6, D and Cq. The 
phase diagram in Fig. 3 is obtained for values of the parameters: Cq = 0.5, a = 0.5, b = 0.4, 
and DVpvK^^ — 20. The membrane shear elasticity is k ~ 3 x lO^^dynecm"^ [0, while the 
plasma density is the same of water: p ~ Igcm"'^; for a flow velocity of the order of a centimeter 
per second, this corresponds to a value I? ~ 6 in Eqn. (7). The drift C(r, z) has been modelled 
putting: C(r, z) = Co(r^ + z^)^^ . It is important to stress, however, that the precise form of C 
is unimportant, provided it does not change sign, which, if occurring in correspondence of the line 
of turning points for 0, could lead to spurious flxed points at (0,a;) ^ (§ + riTr, 0). In the present 
case (see Fig. 3), the distance of the line of turning points from the walls, however, is sufficient to 
guarantee that the far field expression for C(0, x) has at least qualitative value, and that changes 
of sign in C should not be a problem. 

In the situation considered, one finds the transition to a Fahraeus-Linqdwist effect dominated 
situation for d < 5i?. The particular value of d, however is not interesting in itself, depending 
on the choice of the model parameters. What is interesting, is the sharpness of the transition to 
the Fahraeus-Lindqwist dominated regime, which is a consequence of the strong dependence of the 
drift on the parameter r [see Eqn. (8)]. A red cell starting at {9,x) = (0,0.7), while drifting to 
X = r, will travel a distance L along the fiow, that jumps from L ~ 30c?, for r — .25, to L ~ 1500d 
for r — .15. 

To summarize, a qualitative picture arises, of a fiow in a small blood vessel, separated into an 
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Figure 3: Phase diagram of the trajectories in the — x plane for r = 0.2. The dashed Une indicates 

the curve 6 = {] of turning points for 9{t). The trajcc;tory pattern iUustrated depends solely on 
where the line of turning points crosses the x-axis and on where the trajectories bend towards or 
away from the wall, which could be modified only if C presented sign changes. 

inner and outer region, dominated respectively by the quadratic and the linear part of the shear. 
An isolated red cell in the external region, will remain there, in a state of flipping motion for a 
very long time. Once in the internal region, instead, it will dispose itself perpendicular to the flow, 
assuming, under the effect of the quadratic shear, a parachute kind of shape; at the same time, it 
will be pushed towards the vessel axis by interaction with the walls. A Fahraeus-Lindqwist effect 
will therefore occur, if enough of the vessel interior lies in this quadratic shear dominated region. 

In the present study, only the behavior of an isolated red cell has been taken into exam. 
Consideration of the mutual interaction of red cells in real blood, however, is going to strengthen 
the results obtained. Inter-particle interactions will act like a noise in Eqn. (8) and are going to 
produce two effects on the dynamics of a single red cell: the destabilization of the fixed points 
at {9,x) = (f + riTT, 0) and the transverse diffusion of trajectories in the flipping motion region. 
If the quadratic region is wide enough, say, x ~ r^VpuK,~^ > r, the second effect will dominate, 
speeding up collapse of the trajectories onto the channel axis. Thus, a larger values for the value of 
d giving the transition to the Fahraeus-Lindqwist effect dominated regime, is going to be expected, 

with a dependence on the flow parameters given by; d < (Vpv/K^ R, in place of the analogous 
expression provided by Eqn. (9). 
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